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Abstract 

We show that correlators of some local operators in gauge theories are sen- 
sitive to the presence of the instantons even at high temperature where the lat- 
ter are bound into instanton-anti-instanton "molecules". We calculate correlation 
functions of such operators in the deconfmed phase of the 2+1 dimensional Georgi- 
Glashow model and discuss analogous quantities in the chirally symmetric phase of 
QCD. We clarify the mechanism by which the instanton-anti-instanton molecules 
contribute to the anomaly of axial U{\) at high temperature. 
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1 Introduction. 



In this paper we examine one particular aspect of high temperature phase of gauge theo- 
ries. Specifically, many gauge theories at zero temperature have instantons. The behavior 
of the instanton ensemble is crucial for understanding of many dynamical properties of 
these theories. 

In QCD with massless fermions there are indications that the instanton "liquid" is 
responsible for the chiral symmetry breaking [0J. A successful phenomenology of chiral 
symmetry breaking has indeed been developed based on this idea. Another example of 
a gauge theory that supports instantons is the 2+1 dimensional Georgi-Glashow (GG) 
model. In this theory (as well as in 2+1 D compact electrodynamics) the instanton- 
monopole gas has been shown analytically to be responsible for confinement a long time 
ago §. 

When heated above a certain critical temperature these gauge theories undergo a 
phase transition. In QCD the chiral symmetry is restored, while in the GG model the 
deconfining transition occurs. The status of the instantons in the hot phase is less 
certain. One thing is clear - they become less important at high temperature. The 
temperature acts as an infrared cutoff on the instanton size in QCD, suppressing the 
instantons of the size greater than the inverse electric Debye screening mass [[3]. There 
is also another, perhaps more significant effect. In the GG model it has been shown 
analytically |4], || that the instantons at high temperatures are bound into "molecules" . 
The same is believed to be true in QCD 

This "disappearance" of instantons at high temperature has been discussed in the 

4 It has been even suggested that the binding of the instantons drives the chirally restoring phase 
transition This last point has not been definitively proven however. Indeed, this is known not to be 
the case in the GG model |5| , and a physically different mechanism for the chiral transition was suggested 
in Our interest here however is not in the phase transition itself, but rather in the properties of 
the hot phase. In this respect both approaches of Q and agree that the instantons are bound in 
"molecules" . 
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literature. Naively one may expect that once the instantons are bound in pairs, their 
effect in the infrared is negligible. If this is the case, the anomaly of the axial Ua(1) 
symmetry of QCD disappears, and the chiral symmetry of QCD in the hot phase is 
enhanced to SU(Nf) x SU(Nf) x Ua(1). One dynamical effect of such symmetry en- 
hancement would be the degeneracy of the correlation lengths in the scalar (5 and a) 
and the pseudo-scalar (n and rf) channels. In fact various references to the restoration 
of the axial {7a(1) symmetry can be found in the literature |§. 

Of course this argumentation is too naive, and is recognized as such by many authors 
in 0. The statement about the anomaly of Ua(1) in QCD, c^J^ = j^FF, is an 
operatorial one. As an operatorial statement it remains true at nonzero temperature. 
Thus the effects of the anomalous breaking are bound to be seen at some level even at 
high temperature. Indeed the lattice data does show that the scalar and pseudo-scalar 
correlation lengths are not degenerate 0. Some analytic discussion of the consequences 
of the anomaly at high temperature has been given in | ID| . In particular it has been 



shown that in the two flavor case the instantons contribute directly to the correlators 
of the fermionic bilinears, and thus very likely lead to splitting between the scalars and 
the pseudo-scalars. For Nf > 2 no such direct contributions exist. Thus the term 
u Ua(1) restoration" is not to be taken literally. Still since no analytic calculation of the 
instanton effects at high temperature are available, one may have lingering doubts about 
their importance. It would thus be useful to have an explicit analytic calculation of some 
quantities which are most sensitive to the anomaly, and thus exhibit it in a clear way at 
high temperature. 

In QCD this is a difficult task due to the large number of degrees of freedom associated 
with an instanton (color orientations). However in the Georgi-Glashow model situation is 
much simpler, since the instantons are not degenerate. The theory is weakly interacting 
at high temperature, and thus under analytic control. Moreover, in many respects the 
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Georgi-Glashow model is very similar to QCD, especially as far as the role of instantons 
is concerned. The purpose of this paper is therefore to present explicit calculations of 
the instanton effects at high temperature in the Georgi-Glashow model and the close 
analogy between these results and our expectations in QCD. 

The general pattern of symmetries in the theories we have in mind is the following. 
On the classical level the action is invariant under a global symmetry group G. Quantum 
effects break this symmetry anomalously down to its subgroup H. At zero temperature 
H is spontaneously broken. At high temperature H is restored, but G is still broken 
by the anomaly. Thus a correlator O, which is invariant under H, but not invariant 
under G, is allowed to be non-vanishing in the anomalous theory, but has to vanish if 
the anomaly of G were to disappear. This particular set of correlators is therefore likely 
to be very sensitive directly to the instanton effect. We will consider examples of such 
operators in the following. 

There are many parallels between QCD with massless fermions and the 2+1 GG 
model. In particular the confining physics of the latter is in many respects similar to 
the chiral physics of QCD. These parallels been discussed in detail in 0. The basic 
correspondences are the following. 

• Classical axial SU(Nf) x SU(Nf) x Ua(1) symmetry of QCD «-> Classical magnetic 
Um(X) symmetry of the GG model . 

• Axial anomaly due to instantons in QCD<-> Magnetic anomaly due to monopoles 
in the GG model. 

• Non-anomalous SU(N f ) x SU(N f ) subgroup of SU(N f ) x SU(N f ) x U A (1) in QCD 
<-> Non-anomalous Z 2 subgroup of Z7m(1) m the GG model. 

• Spontaneous breaking of SU(Nf) x SU(Nf) by the chiral condensate < ipip > in 
QCD^ Spontaneous breaking of Z 2 by the vortex condensate < V > in the GG model. 
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Keeping this analogy in mind, we proceed to analyze in the next section the monopole- 
instanton effects in the 2+1D GG model at high temperature. Since the theory is weakly 
interacting, the calculations here are under complete analytic control. In section 3 we 
calculate explicitly multi-instanton contributions to relevant correlation functions. Our 
results illustrate very clearly the effects that we also expect to take place in QCD. We 
conclude in section 4 with discussing further the analogy with QCD. 

2 The instantons and their interactions in the Georgi- 
Glashow model at high temperature. 

2.1 The model. 

Consider the SU(2) gauge theory with a scalar field in the adjoint representation in 2+1 
dimensions. 



±(D,h a ) 2 + ^(h a h a -v 2 ) 2 



(1) 



2<? 

In the weakly coupled regime v ^> g 2 , perturbatively the gauge group is broken to 
U(l) by the large expectation value of the Higgs field. The photon associated with the 
unbroken subgroup is massless whereas the Higgs and the other two gauge bosons W ± 
are heavy with the masses 

M 2 H = 2Xv 2 , M 2 W = g 2 v 2 . (2) 

Thus perturbatively the theory behaves very much like electrodynamics with spin one 
charged matter. 

This theory has a global symmetry - the magnetic symmetry JOJ [T^| . Classically the 
following gauge invariant current is conserved 

F» = Py a - -^ vX e abc n a {V u n)\V x n) c (3) 



with n a a unit vector in the color space in the direction of the Higgs field, n a = h a /\h\. 
This defines a conserved charge through $ = J d 2 xB. This continuous symmetry is 
spontaneously broken in the vacuum, and the massless photon is the Goldstone boson 
which reflects this breaking in the spectrum. 

However as is well known the classical action eq.(|l|) supports instanton-monopole 
solutions with finite Euclidean action 

h a (x) = x a h{r) 

Al{x) = l - [e aixv x v (l - 00 + 5 a y 2 + (rA - 2 )x%] (4) 

where x a = x a jr. In the presence of such a monopole the magnetic current is not con- 
served, but rather has a non- vanishing divergence proportional to the monopole density. 

d»F, = jp (5) 

The Um{^) magnetic symmetry is thus anomalous in the quantum theory. As discussed 
in detail in |12| only the discrete Z 2 subgroup is unaffected by anomaly and thus remains 
a symmetry in the full quantum theory. 

The local order parameter that transforms as an eigen-operator under the Z 2 magnetic 
symmetry is the magnetic vortex operator V: 

V -> -V (6) 



The operator V has been constructed explicitly in JIT] , [L2 

V(x) = exp - / d 2 y 
a J 



9 

= exp — / dy % e ijn a E*(y) (7) 
9 Jc J 

One can think of it as a singular 577(2) gauge transformation with the field dependent 
gauge function 

X a (y) = -&(x - y)n a (y) (8) 
9 
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The dynamical effects of the anomaly are twofold. First, the photon becomes a 
pseudo-Goldstone boson and acquires a finite albeit small mass m 2 h oc exp{— At[M w / g 2 } 
(in the BPS limit) . Another effect is confinement of bosons with the string tension 
a oc g 2 m ph . 

2 

When heated the model undergoes a phase transition to deconfined phase at T c = |- 
H . At this transition the non-anomalous magnetic Z 2 symmetry is restored. 

2.2 The effective action and the instanton interactions. 

We are interested primarily in the high temperature phase, T > T c . At these tempera- 
tures the dimensional reduction applies and the effective Lagrangian was derived in |5j 
in terms of the abelian Polyakov line P. 




Here the Polyakov line P is a pure phase P = exp{i^Ao}, and Aq is the component of 
the vector potential parallel to the direction of the Higgs field n|. 

The action eq.(^) deserves some explanation. Firstly, as any dimensionally reduced 
Lagrangian it is valid at distance scales greater than the inverse temperature. The 
parameter a appearing in it is the ultraviolet cutoff, and is of this order a = 2%/T. 
The second term expanded to order Aq gives the electric Debye screening mass to Aq. 
The mass squared is proportional to the fugacity of the charged particles, which in this 
model are the bosons Mf, cx exp{— Mw/T}. The exact proportionality coefficient 
is straightforwardly calculable at one loop, but is inessential for our purposes. The 

last term is induced by the monopole-instantons. As discussed in || the monopoles 

5 Note that the transition temperature is much less than the Higgs expectation value T c « v 2 . Thus 
the direction of the Higgs field is well defined also above the transition temperature. Of course at high 
enough temperature T = 0(v 2 ) the Higgs field will also fluctuate strongly and the nonabelian part 
of the gauge group will also become important. At these temperatures the theory becomes essentially 
nonabelian with all the ensuing complications. To keep our calculation under analytic control we restrict 
ourselves to temperatures T c <T < v 2 . 
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in this effective action appear as vortices of the Polyakov line P with unit winding. 
Indeed the last term in eq. (Q) is of the type of Skyrme term and vanishes on any smooth 
configuration of A$. Its only contribution comes from points at which the winding of 
P is concentrated. It therefore just counts the number of vortices and anti-vortices of 
P, that is of monopoles and anti-monopoles. On a configuration with one monopole the 
contribution of this term to the partition function is equal to (. Thus ( is the fugacity of 
the monopole (or anti-monopole), ( oc exp{— 47rM w /g 2 } []. An important point is that 
this dimensionally reduced action has not been derived in the derivative expansion, but 
rather in the expansion in powers of the exponentially small parameters M^a and C, H . 
It is thus valid for all distance scales above the UV cutoff a. 

It is easy to see from the effective action eq.@ that the monopoles are bound by a 
(screened) linear potential. Consider a configuration with unit winding of the field P. 
Due to the potential term P 2 the minimal action configuration can not be a rotationally 
symmetric hedgehog. Such a configuration would "cost" action proportional to the vol- 
ume, since the field P would be away from its vacuum value everywhere in space. The 
best one can do is to have a quasi one dimensional strip in which the winding is concen- 
trated, while everywhere else in space P would be equal to 1 (or —1). This configuration 
is schematically depicted on Fig. 1. The width of the "dual confining string" must be 
clearly of order of Mp , while the action per unit length s oc p-M^. Thus the action of a 
single monopole diverges linearly with the size of the system. Obviously the monopole- 
anti- monopole pair separated by a distance L have the action Ls. When the distance 
is large enough, another pair can be popped from the thermal ensemble and screen the 
linear potential. The critical distance at which this happens is clearly determined byQ 

8ttM w 

L c s = 21ogC = y~ 10 
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6 The proportionality constant is slightly different than at zero temperature reflecting the fact that 
the field modes with momenta p > T have been integrated out to arrive at the effective action eq.(^)[Q. 

7 In this equation we neglected the effects of the finite temperature "renormalization" of the monopole 
fugacity C- 
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Figure 1: The string- like configuration of the Polyakov line P that accompanies the 
monopole-instanton in the high temperature phase. 



or 

Lc^Mj (11) 

Thus as long as the temperature is much lower than M^, the length of the "dual string" 
is much greater than its thickness. The "potential" between the monopoles is therefore 
linear, but screened at large distances, much in the same way as the confining potential 
in gauge theories with heavy fundamental charges^. 

Now it is interesting to consider the correlation functions of the vortex operators. 
There is a direct analogy between these correlation functions and the correlators of the 
fermionic bilinears in QCD. In particular if the magnetic U^(l) were restored at high 
temperature, this would imply a simple relation between correlation lengths in different 

channels, directly analogous to the would be degeneracy between the rj and the pions in 

8 Of course we should always keep in mind that while talking about instanton "potential" we really 
mean "action" and not energy. Thus physically "confinement" of instantons is very different that 
confinement of charged particles. 
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QCD. Recall that the vortex operator in the GG model is a pseudo-scalar [|12j, so that 
the parity transformation acts on it as 

V^V* (12) 

The parity even and odd eigenstates are therefore 

V + V* V - V* 

V + =^, V~=i— (13) 

If the magnetic £/m(1) is restored, the correlators < VV > and < V*V* > both vanish, 
and thus the following relation should hold 

< V + (x)V + (y) >=< V-(x)V~(y) > (14) 

The correlation lengths in the parity even and the parity odd channels should be the 
same. This is the direct analog of the would be equality of the correlation lengths in 
the 7] = tpjsip and the a = ipip channels in QCD. We will see in the following that this 
equality in fact does not hold due to the instanton effects. 

3 Instanton effects at high temperature. 
3.1 Um(X) invariant correlators of V. 

Let us remind the reader how calculation of the correlators of V proceeds in the frame- 
work of the effective action eq.(§). For pure gauge theories this has been discussed in 



detail in fLo| . In the present, essentially Abelian case the procedure is very similar. The 
vortex operator V(x) in eq.(^) induces a singular abelian gauge transformation which 
forces the vector potential A to jump by n across the curve C in the definition eq.([7|). 
Thus the Polyakov line P is forced to change sign across C. Formally, the insertion of 
the operator V(x) into the Euclidean finite temperature path integral leads to the shift 

diA -> diA - cti (15) 
10 




Figure 2: The configuration of the Polyakov line P that dominates the correlator 
V*(x)V(y). 



in the derivative terms in eq.(H). Here a, is a "vector potential" of a point-like vortex of 
vorticity 1/2 at the point x 



a i (y)=7m^(y)5(y-C) (16) 

where n (y) is a unit vector normal to the curve C at the point y. The result does not 
depend on the curve C, but only on its end point x, since the shape of the curve may be 
changed into C by transforming P — > —P in the area enclosed by C — C'\IE\. 



The calculation of the correlation function V(x)V*(y) leads to the path integral 
with the vortex at x and the anti-vortex at y. At high enough temperature due to the 
factor T/g 2 in the action this path integral can be calculated in the steepest descent 
approximation. The integral then is dominated by the solution of the classical equations 
of motion subject to the condition that P has to change sign across the curve C that 
connects the points x and y, see Fig. 2. 



The shape of this solution is easy to understand [15|. In the hot phase the configu- 
rations P = 1 and P = — 1 are degenerate. There are therefore classical solutions which 
have the form of domain walls, that interpolate between these two configurations, the 
so called Z 2 domain walls [|T6|]. Since our effective action is a simple sine-Gordon theory, 
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the solution for such a domain wall is known explicitly^ 



1nT 

A = arctane^"^ (17) 
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The action of this solution per unit length is 



~ 8T 
cr= — M D 

g 2 



In the calculation of V(x)V*(y) the discontinuity in the boundary conditions forces 
P to jump from the vicinity of one "vacuum" to the vicinity of the other across C . But 
since P has to go to the same value in all directions at infinity, it has to return to the 
original vacuum. The solution thus is a domain wall of Fig.2. Using the result eq.(|18D 
we find that for large separations^ 

< V(x)V*(y) >oc exp{-|x - y\a} (19) 

Actually one can do a little better and determine the prefactor quite easily. The prefactor 
originates from the action associated with the endpoints x and y at which the wall 
terminates. Close to the endpoints (at distance smaller than M^ 1 ) the mass term is 
irrelevant and the configuration of Aq must be that of a rotationally symmetric vortex 
(anti-vortex) of vorticity 7r 

M*) = (2°) 

with 9 a planar angle relative to the location of the vortex. Such a vortex carries the 
Coulomb "self energy" 



T c^d 1 dr ttT 

2f ^ la 7" ~ ~ V 



5S = —2^I — = -— MaM D ) (21) 



9 The Skyrme term in eq.(^|) does not play any role in this discussion, since the solution in question 
is smooth and this term vanishes for smooth functions Aq. 

10 Again the Skyrme term does not contribute, since the contribution of the cut C is subtracted 
explicitly by the shift eq.(|l5|). 
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Adding the contribution of both end points we fincQ 

< V{x)V*(y) >= (aM D )^ exp{-\x - y\a} (22) 

Note that this calculation was performed entirely in the zero instanton sector. There 
is of course an instanton-anti-instanton contribution. We will discuss it at the end of 
this section. 

3.2 Instanton sensitive correlators. 

Let us now consider the calculation of < V 2 (0) >, which is the simplest operator not 
invariant under the magnetic C/jw(1), but still invariant under the non-anomalous Z2. 
Neglecting the instanton contributions, we should find that the VEV of this operator 
vanishes. Indeed the path integral for < V 2 (0) > has the boundary condition imposing 
vorticity 2tt of P at the point x = 0. In the zero instanton sector the action with 
such boundary condition diverges. The reason is precisely the same as the one for linear 
potential between the instantons. The external vortex associated with V 2 (0) will pick up 
a "dual confining string", which carries finite action per unit length. This configuration 
has action proportional to the linear size of the system and thus the VEV of V 2 vanishes 
exponentially in the thermodynamic limit. The leading non- vanishing contribution to 
< V 2 (0) > comes from the sector with one anti- instanton . Clearly if there is an anti- 
instanton in the vicinity of the point x — 0, the dual string that originates at x = 
will terminate on it and thus the result will be finite. To calculate the expectation value 
we note that it is dominated by the anti-instanton sitting relatively close to the point 
x = 0, since otherwise the action of the configuration is very large. On the other hand if 

the coordinate of the anti-instanton x a is smaller than M^ 1 , the action of the classical 

11 In principle we should also subtract 2Af£ ) 1 from the length of the string for consistency. This 
correction is however sub-leading, since it brings in a factor exp{2aM^ 1 } which is of order 1 in the 
interesting temperature range. 
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solution is just the action of a dipole with the vector potential 



ao(x) = 9{x) — 6{x — x a ) 



(23) 



The expectation value is therefore given by the integral over the quasi zero mode x a 



Let us now consider the correlation function < V(x)V(y) >. As far as the magnetic 
quantum numbers are concerned, this correlator is the same as the composite operator 
V 2 , but it is more interesting since it must also exhibit some physical correlation length. 
The calculation of this quantity is not much different. In the zero instanton sector it 
vanishes for the same reason as V 2 . It carries a net "vorticity" and this vorticity has 
to be screened by an anti-instanton in order to give a finite result. In principle the 
anti-instanton can be located at an arbitrary point in space. The field configuration 
corresponding to two external vortices associated with insertions of V and the anti- 
instanton is schematically depicted on Fig.3. 

Clearly the contribution of anti-instantons sitting at a significant distance from the 
straight line connecting x and y is exponentially suppressed, since the "dual confining 
string" in this configuration is longer. Thus the leading contribution is given by the anti- 
instantons sitting within the transverse distance M^ 1 of the straight line. The correlation 
function thus is 



where we have taken x = (xi,0) and y = (yi,0). Here z is the position of the anti- 
instanton. As long as \x — y\ » M^ 1 , the main contribution comes from the anti- 
instantons far from the endpoints x and y. For these configurations the integral over z\ 




(24) 



< V(x)V(y) >= C / -4e W {S(x,y,z)} = ( f 

J 7TGT Jx 



y dzi 



nM D a 2 



exp{-S(x,y, Zl )} (25) 
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Figure 3: The configuration of one anti-instanton contributing to the correlator 
V(x)V(y). 



just gives \x — y\, and the action is the sum of the action of the string a\x — y\ plus the 
self energies of the anti-instanton and the vortices at the endpoints x and y. Those are 



calculated as in eqs.( pT| , |22D . So that finally we get 



< V(x)V(y) >= C{aM D )^ — — - exp{-|z - y\a} (26) 

d 

Interestingly enough the correlation length that appears in this correlation function is 
the same as the one appearing in the correlator of vortex and anti- vortex V(x)V*(y). 

Thus as expected, both the VEV of V 2 and the correlator V(x)V(y) are sensitive to 
the instantons. Both these quantities are indeed directly proportional to the fugacity 
of the instanton (. Even though in the £7m(1) invariant sector the instantons appear 
only in bound pairs, in the calculation of the £7m(1) non-invariant correlation functions 
they do appear separately. Such a single instanton (anti-instanton) binds to the external 
insertion that carries the explicit Um(1) charge and screens this charge. This leads to 
non- vanishing correlation functions in the channels whose quantum numbers match those 
of an integer number of instantons (anti-instantons). 
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3.3 Multi-instanton contributions. 



As mentioned earlier, there are also multi-instanton contributions to the vortex corre- 
lation functions. Those become important when the separation \x — y\ is large enough. 
The reason is that although each extra instanton-anti-instanton pair is suppressed by 
( 2 , at large \x — y\ the "entropy" is large enough and thus these contributions are not 
negligible. 

Consider first the correlator V(x)V*(y). In order that the multi-instanton contri- 
bution not be suppressed exponentially, all instantons must be located on the straight 
line connecting the points x and y. Moreover, the instantons and anti-instantons must 
alternate, since otherwise some sections of the dual confining string connect regions of 
vorticity 3, and thus has much higher dual string tension. Other than these constraints, 
the coordinates of instantons and anti-instantons are arbitrary and have to be integrated 
over. The contribution of n instantons and n anti-instantons is therefore given by 

TvnT 

<V(x)V*(y)> n = (aM D )~P- exp{-\x - y\a} (27) 

i-xnT 9 -2n f V f Vn f Xn f V1 

x C n (aM D )^7cM D a 2 / dy n / dx n / dy n -x... / dx x 

Here Xi and are longitudinal coordinates of instantons and anti-instantons correspond- 
ingly 0. Performing the integration and summing over n we find 

< V(x)V*(y) >= {aM D )^ exp{-|a; - y\d} cos\i{({aM D )l^~ l ^ X - ^ } (28) 

The multi-instantons have therefore a very nontrivial effect on the vortex correlator. 
They "split" the correlation length in two, generating two distinct correlation lengths 

a± = a ± C{aM D )^~ l — (29) 

7ia 



12 Note that the usual symmetry factor does not appear in this expression due to ordering of the 
coordinates of the instantons and anti-instantons. 
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The calculation of the multi-instanton contributions to V(x)V(y) is very similar, and 
gives 

ttT 2tiiT -i \X — 7/ 

< V(x)V(y) >= (aM D )^ exp{-|x - y\a} smh{((aM D )^ ± ^} (30) 

Again we see that the same set of correlation lengths appears in the two correlation 
functions. 

Eqs.( |28| , |30"l) are very simple from the point of view of the parity eigenstates V + and 



V eg. ([13]) . They can be rewritten as 

< V + (x)V + (y) > = (aM D )^~ exp{-\x - y\a-} 
<V-(x)V~{y)> = (aM D )^ exp{-|x - y\a + } (31) 

Thus we see explicitly that the correlation lengths in the scalar and pseudo-scalar chan- 
nels are not equal due to the instanton effects. This difference decreases with temper- 
ature, but in the temperature range under consideration (T <C My/) always remains 
finite. Recalling that cx exp{— Mw/T} and ( oc exp{— 4nMy//g 2 } we find that the 
ratio of the difference of the inverse correlation lengths to their sum varies between 

r vM w g 2 

~ , - = exp{ — }, T^T C = — (32) 

er+ + cr_ g z Air 

and 

- — — — = exp{ — }, T -> M w (33) 

o-+ + a- g 2 

Finally it is worth mentioning that the quantities that get the direct instantons con- 
tributions naturally are proportional to the instanton fugacity, but their smallness is not 
uniform. For example as already mentioned, the correlation lengths in the vortex-vortex 
channel is the same as in the vortex-anti-vortex channel. Moreover comparing the vortex 
- vortex correlator eq.(ftDD with the vortex-anti- vortex correlator eq.(P%D we see that al- 
though the former is smaller at short distances , at distances of order 7ra£ -1 (aMD) 
the two are numerically practically equal. 

17 



2nT 

g 2 



3.4 N > 2. 



Let us briefly comment on the situation in the GG type models for N > 2. The quan- 



titative discussion of the phase transition will be given elsewhere p7| . Here we want to 
discuss only the main qualitative features. The basic degree of freedom is still a vortex 
operator V(x), and the effective theory at high temperature is similar to eq.(^) with the 
main difference that the potential term is P N + P* N , and thus is invariant^. The 
instanton again is just a vortex of P with unit vorticity. 

First consider the calculation of < VV* >. Just like for N = 2, the contribution in 
zero instanton sector is given by the Zn domain wall 

< V(x)V*(y) >oc exp{-|x - y\a} (34) 

The instanton-anti-instanton contribution now is however different. The point is that 
the instanton (anti- instanton), having winding number one, is the source (sink) of N 
Zn domain walls. Thus the system of the vortex V(x) plus an anti-instanton at some 
point X\ together are the sink of N — 1 domain walls. The string tension of this bunch 
is roughly (N — 1)5". Thus in order for the instanton-anti-instanton contribution not to 
be suppressed by the exponential factor, the anti-instanton must sit within a distance 
[(N - 2)<t]- 1 of the instanton, see Fig. 4. 

Again summing over the number of the AI pairs, we find that the contribution expo- 
nentiates 



< V(x)V*(y) >= (aMn^expi-lx-ylajexpieiaMn)^- 1 ' y] — } (35) 

7ra z {T\ — l)o 

13 The actual situation generically is more complicated than that. In particular there are N— 1 distinct 
vortex operators. Also there are N — 1 distinct Abelian Polyakov lines, corresponding to the N — 1 
Cartan algebra generators, and all of those should appear in the effective action. However one can 
always choose the parameters of the model to be such that one of these operators is much lighter than 
the others. The rest of the Polyakov lines then can be integrated out and the effective action depends 
only on one P. For simplicity we consider only this situation in the present article, although the main 
conclusions are the same also in the generic case. 



18 



t ♦ t t t I V 1 1 1 1 1 V f f ) H f f ♦ I 



x 



x A 



x 



x. x. x. x^ 



X X X X 
x x X x 



X 



■t 



t t i 



X X x X 



XXX 



(Ml 



AAAAAAAAAAAAAAAAAAAAAAAAA 

Figure 4: The instanton-anti-instanton configuration contributing to the correlator 
V(x)V*(y) for A = 3. The instanton (anti-instanton) is a source (sink) of three ele- 
mentary Z3 domain walls. The curve C is the "fake" infinitely thin wall as in eq. (|I~6D . 



So the result is quite distinct from the A = 2 case. Only one correlation length is 
present, although there is indeed a direct contribution of instantons into this length of 
order 0(C 2 )- 

The behavior of the correlator < V(x)V(y) > is also very different. Since this corre- 
lation function is not symmetric under the magnetic symmetry, which is unbroken 
in the high temperature phase, it strictly vanishes. Instead, the instantons contribute 
to the correlation functions with the quantum numbers of V . In particular consider 

< V(x)V N ~ 1 (y) >. The leading contribution to this correlation function comes from the 
configuration with an anti-instanton within the distance [(A — 2)a]~ 1 of the point y and 
an arbitrary number of the AI pairs. 

< V{x)V N -\y) >= aaMn)^^- 1 exp{-|x - y\a} exp^M^' 1 J^L }(36) 

7rcr(A — 2)a 

In this calculation we only kept the leading contribution to the correlation function, 
which in the absence of instantons decays with the correlation length l^ 1 = 5. There 
is of course another contribution. The configurations where the instantons and anti- 
instantons are allowed to separate will lead to the component with the correlation length 
/ _1 = (A — 1)5". This is precisely the correlation length that in the absence of instantons 
dominates the correlator < V N ~ l {x)V* <yN ^ 1 \y) >. The instantons clearly lead to the 
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mixing of the operators V and V*^ 1 - ) in full analogy with the case N = 2. This mixing 
is accompanied with the instanton induced increase in the longest correlation length 



-i 



[a - (aM D ) * 



-1 



c 2 



-1 



(37) 



wa 2 (N-2)a 



Thus for N > 2 the instantons do not remove the degeneracy between V + V* and 
V — V*,but they do significantly affect other correlation functions consistently with the 
magnetic symmetry. 

We conclude that even though at high temperature the instantons are bound in pairs, 
they are still relevant in the infrared. The anomaly does not disappear. Its manifestation 
is non-vanishing of non-singlet correlators as well as the difference of the correlation 
lengths in the scalar and pseudo-scalar channels for N = 2 and the operator mixing 
consistent with the unbroken magnetic for N > 2. 

4 Discussion 

The calculation presented in this paper illustrates two simple points. First, which cor- 
relation functions vanish and which don't is determined entirely by the non-anomalous 
symmetries. Second, the instanton contributions are important at high temperature even 
though the instantons in the partition function are bound in pairs. Both these features 
we expect of course to hold also in QCD. 

In particular let us consider the fermionic bilinear correlation functions. The fermionic 
bilinears that transform as irreducible representations of the non-anomalous SU(Nf) x 
SU(Nf) are 



(38) 
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whereas the parity eigen-operators are 

0± = 0±0 ] (39) 

Now consider the correlation functions + (x)0 + (y) and 0-(x)0-(y). If the corre- 
lators T =< 0(x)0(y) > and < T* = 0*(x)0*(y) > both vanish, then the correlation 
functions of the parity even and the parity odd operators are equal. Whether this is 
true or not depends on the number of flavors. For N > 2 the vanishing of these corre- 
lation functions is assured by the SU(Nf) x SU(Nf) symmetry, since T is the product 
of two fundamental representations (separately for the right and left symmetries) and 
therefore does not contain a singlet. For Nf = 2 on the other hand the product of two 
fundamentals does contain a singlet. Thus there is no symmetry obstruction for T to be 
nonzero^. 

Dynamically the contributions to T come from the instanton sector. For Nf = 2 the 
number of fermionic operators in T is just right to saturate the fermionic zero modes 
of an instanton. Recall that the fermionic zero modes at high temperature away from 
the core of the instanton exponentially decrease as exp{— 7rT|x|} |IJ. This is the 
same exponential decay as of the lowest fermionic Matsubara mode. Thus we expect 
the situation to be very similar to the N = 2 Georgi-Glashow model. The perturbative 
behavior of the correlator 0(x)0*(y) is 

< 0(x)0\y) >oc exp{-27rT|x - y\} (40) 

The instanton-anti-instanton contributions have the same exponential behavior, and 

extra powers of \x — y\ due to the arbitrariness of the instanton's (anti-instanton's) 

14 In this discussion we did not mention at all the anomalous J7a(1). Of course, if we assume that it 
is restored, the parity odd and even channels must be degenerate, since they transform into each other 
under the axial U(l). 
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position. It is very likely that the sum over the number of instantons exponentiates as 
in eq. (ffif ) and results in the splitting of the basic correlation length into two 2nT ±x(. 
These two distinct correlation length then appear separately in the parity even and parity 
odd channels. 

For Nf > 2 the number of fermions in T does not match the number of the fermionic 
zero modes in the instantons, and thus the only way instantons can contribute is through 
AI pairs. They will however give direct contributions to the correlators of the type 
T N =< 0(x)0 2 N-2(y) > where 2 n-2 is the operator containing 2N — 2 fermions, such 
that the quantum numbers of T/v are the same as of 't Hooft's determinant. This is 



essentially the argumentation of |ID |. Again this is mimicked perfectly by our discussion 
of the N > 2 GG model. 

Although analytic multi-instanton calculations in QCD are difficult, the toy GG 
model gives us a very clear picture of what the expected result is. Hopefully this knowl- 
edge may give us a hint of how to separate the most important contributions so that an 
analytic calculation may become possible after all. 
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